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computational  grid  by  generating  aggregates  of  rigid,  Immobile,  and 
Impermeable  flow-field  cells.  The  simple  transmissive  boundaries  in  HULL 
simulate  a  zero-gradient  condition  across  the  boundary  for  both  Che  pressure 
and  Che  normal  component  of  velocity.  Simple  transmissive  boundaries  such  as 
these  will  partially  reflect  waves  that  strike  them.  Including  shock, 
compression,  and  expansion  waves.  The  strength  of  these  reflected  waves  Is 
directly  related  Co  Che  strength  of  Che  Incident  waves.  These  reflected  waves 
Chen  travel  back  Into  Che  computational  grid,  modifying  the  flow-field 
conditions  In  Che  regions  through  which  they  pass,  thereby  ending  Che 
simulation  of  free-field  condltlons.^The  usual  procedure  is  to  locate  the 
transmissive  boundaries  far  enough  amy  from  the  target  so  that  these 
artificial  waves  do  not  arrive  at  the  target  during  the  time  of  interest. 
However,  limitations  due  to  available  computer  storage,  cost,  speed,  or  any 
combination  thereof,  may  necessitate  a  computational  grid  design  such  chat 
some  disturbances  at  Che  target  frod  Che  transmissive  boundaries  are  probable. 
Hence  It  Is  necessary  to  estimate  the  effects  of  these  disturbances  on  the 
computed  loading  of  Che  target.  The  analysis  of  Che  computations  performed 
here  provides  Insight  Into  the  nature  of  the  disturbances  and  a  means  for 
estimating  Che  deviation  from  free-field  results. 
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I.  INTRODUCTION 


One  of  the  principal  uses  of  the  US  Army  Ballistic  Research  Laboratory 
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(BRL)  version  of  the  HULL  *  hydrocode  has  been  to  perform  blast  loading 
computations  on  military  targets.  The  three-dimensional  (3-0)  version  of  HULL 
has  been  shown  to  give  excellent  agreement  with  experimental  results  for 
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diffraction  loading  on  a  box-like  target  for  a  weak  shock  at  obliquity  . 
However,  that  study  was  carried  only  through  the  shock  diffraction  phase. 

4 

More  recent  studies  have  been  carried  out  using  the  cyllndrically  symmetric 
version  of  the  HULL  hydrocode,  with  variation  axially  and  radially  but  not 
rotatlonally.  The  target  used  in  chose  studies  was  a  rigid,  right-circular 
cylinder,  located  coaxially  within  a  computational  flow  field  simulating  a 
cylindrical  shock  Cube.  Computations  were  carried  out  for  non-decaying  (step) 
shocks  of  three  different  strengths  for  several  shock  crossing  times  so  as  to 
estimate  Che  effect  of  shock  tube  blockage  on  drag  phase  loading.  In  these 
computations,  Che  outer-radial  boundaries  of  the  shock  tube  were  perfectly- 
reflecting  in  order  to  simulate  Che  interaction  of  wave  reflections  between 
Che  target  and  shock  tube  walls.  This  work  was  ended  in  another  computational 

study^  and  a  combined  computational  and  experimental  study^. 

The  conditions  for  a  reflective  boundary  are  reasonably  well-defined  by 
creating  mirror  images  in  virtual  cells  beyond  Che  active  flow  field.  The 
pressure  in  Che  virtual  cell  is  set  equal  to  that  of  the  adjacent  flow-field 
cell,  and  the  normal  component  of  velocity  is  set  equal  in  magnitude  and 
opposite  in  sign.  The  simple  transmissive  boundaries  in  HULL  simulate  a 
zero-gradient  condition  across  the  boundary  for  both  the  pressure  and  the 
normal  component  of  velocity.  Simple  transmissive  boundaries  such  as  these 
will  partially  reflect  waves  chat  strike  them,  including  shock,  compression, 
and  expansion  waves.  These  waves  typically  reflect  in  kind  from  these 
transmissive  boundaries,  with  their  strength  directly  related  to  the  strength 
of  Che  incident  waves.  These  reflected  waves  Chen  travel  back  into  Che 
computational  grid,  modifying  Che  flow-field  conditions  in  the  regions  through 
which  they  pass,  thereby  ending  Che  simulation  of  free-fleld  conditions.  The 
usual  procedure  is  Co  locate  the  transmissive  boundaries  far  enough  away  from 
the  target  so  Chat  these  artificial  waves  do  not  arrive  at  Che  target  during 
the  time  of  interest.  However,  limitations  due  to  available  computer  storage, 
cost,  speed,  or  any  combination  thereof,  may  necessitate  a  computational  grid 
design  such  that  some  disturbances  at  Che  target  from  Che  transmissive 
boundaries  are  probable.  This  report  presents  the  results  from  a  systematic 
study  of  boundary  effects  for  a  2-D  target  of  modest  complexity  struck  by  a 
strong,  non-decaying  shock  wave. 


II.  TARGET 

The  nominal  target  used  in  this  study,  shown  in  Figure  1,  is 
two-dimensional  and  has  a  "barn-like"  cross-section.  It  has  a  height  of  1.83 
m  (6.0  ft)  and  a  depth  (in  the  direction  of  shock  travel)  of  3.05  m  (10.0  ft). 
The  Cop  face  of  the  target  has  front  and  back  sloping  sections,  each  at  an 
angle  of  30  degrees  to  the  horizontal,  and  a  central  flat  section.  The 
vertical  sections  of  the  front  and  back  faces  are  each  1.22  m  (A.O  ft)  high. 
There  are  two  variations  on  this  target.  The  first,  hereinafter  referred  to 


as  target  SS  (for  Stair-Stepped);  uses  only  fully  rigid  (Island)  cells  to 
construct  the  target,  so  the  front  and  back  slopes  for  the  top  face  actually 
resemble  "stair-steps."  The  second,  hereinafter  referred  to  as  target  RT  (for 
Ramped  Target),  uses  combinations  of  Island  cells  and  half-rigid, 
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half-hydrodynamic  (shore)  cells  to  construct  the  target.  Most  of  the 
computations  were  done  using  target  SS  because  the  shore  cells  are  available 
only  in  2-D  HULL  and  they  cause  the  code  to  run  somewhat  slower.  The 
computations  for  target  RT  were  done  to  provide  an  Indication  of  the  effect  of 
using  the  shore  cells  to  simulate  a  smooth  ramp  as  opposed  to  using  a  stepped 
approximation.  A  reference  computation  that  was  hydrodynamically 
two-dimensional  was  performed  for  target  SS  using  the  3-D  HULL  coding.  In 
this  way.  Information  gained  from  the  2-D  study  can  be  linked  to  similar  3-D 
computations.  For  the  remainder  of  this  report,  "front  face"  refers  to  the 
vertical  section  of  the  front  of  either  SS  or  RT,  "back  face"  to  the  vertical 
section  of  the  back,  "front  slope"  to  the  sloping  section  of  the  top  face, 
"back  slope"  to  the  back  sloping  section  of  the  top  face,  and  "top  flat 
section"  to  the  flat  section  in  the  center  of  the  top  face  (see  Figure  1). 


III.  SIMULATED  CONDITIONS 

The  ambient  conditions  were  chosen  to  be  characteristic  of  those  found  at 
the  White  Sands  Missile  Range.  The  ambient  conditions  were: 

static  pressure  =  85.5  kPa  (12. A  psi), 

static  density  =  0.999  kg/m^  (0.0624  Ibm/ft^), 

static  temperature  =  298.15  Kelvins  (536.7  degrees  Ranklne), 

local  sound  speed  =  346.1  m/s  (1136  ft/s). 

The  blast  wave  considered  is  from  a  206.8  kPa  (30  psi)  overpressure  step 
shock.  The  shocked-gas  conditions  were: 

static  pressure  =  292.3  kPa  (42.4  psi), 

static  density  =  2.282  kg/m^  (0.1425  Ibm/ft^), 

static  temperature  =  446.3  Kelvins  (803.4  degrees  Rankine), 

local  sound  speed  =  423.5  m/s  (1390  ft/s), 

local  Mach  number  =  0.806, 

dynamic  pressure  =  132.8  kPa  (.19.26  psi), 

shock  speed  =  606.9  m/s  (1991  ft/s), 

shock  Mach  number  =  1.753, 

gas  velocity  =  341.2  m/s  (1119  ft/s). 

The  theoretical  conditions  on  the  target  front  vertical  wall  for  normal 
reflection  were; 

static  pressure  =  817.9  kPa  (118.6  psi), 
static  overpressure  =  732.4  kPa  (106.2  psi), 
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static  density  =  4.613  kg/ra  (0.2880  Ibm/ft  ), 
static  temperature  =  617.6  Kelvins  (1112  degrees  Rankine), 
local  sound  speed  =»  498.2  m/s  (1634  ft/s), 
shock  speed  675.1  m/s  (2215  ft/s). 


IV.  HYDROCODE  COMPUTATIONS 


A.  General  Information. 

The  computations  were  performed  on  the  BRL's  CDC  7600  using  the  BRL 
1  2 

version  *  of  the  HULL  hydrocode,  a  modified  version  of  one  obtained  from  the 
US  Air  Force  Weapons  Laboratory  in  1978.  The  code  solves  the  inviscid  Euler 
equations,  using  an  explicit  time  step  finite  difference  algorithm 
(essentially  a  modified  Lax-Wendroff  computing  scheme).  Three  series  of 
two-dimensional  computations  were  performed  for  this  study.  The  first  set  was 
performed  as  a  specific  check  on  the  influence  of  artificially  reflected  waves 
from  the  top  transmissive  boundary  on  the  loading  of  target  SS.  The  second 
set  of  computations,  one  of  which  is  discussed  here,  showed  the  Influence  of 
the  positions  of  the  upstream  and  downstream  transmissive  boundaries  on  the 
loading  on  target  SS.  The  final  set  was  done  as  a  reference  set  to  show  the 
effects  of  using  smooth  slopes  (target  RT)  versus  stair-stepped  slopes  (target 
SS)  on  the  computed  loading  of  the  nominal  target.  This  final  set  also 
Included  a  hydrodynamlcally  2-D  computation  using  the  3-D  coding  to  provide  a 
check  on  the  Internal  consistency  within  HULL,  as  the  2-D  and  3-D  coding  have 
some  differences  between  them.  All  times  quoted  in  these  studies  have  a 
reference  value  of  zero  for  the  theoretical  arrival  of  the  incident  shock  at 
the  front  face  of  the  target. 

All  computations  were  performed  with  the  bottom  boundary  of  the  grid 
defined  as  a  perfect  reflector,  with  no  possibility  of  flow  through  it.  The 
upstream  boundary  was  defined  as  a  constant  Inflow  boundary,  with  input  flow 
conditions  equal  to  those  for  the  steady  206.8  kPa  overpressure  shock  defined  in 
I  the  previous  section.  Any  wave,  real  or  artificial,  arriving  at  this  boundary 

does  not  directly  affect  the  fixed  boundary  conditions  there,  but  will 
typically  reflect  some  manner  of  artificial  wave  back  into  the  computational 
grid.  Hence,  any  reference  in  the  following  text  to  artificial  waves 
originating  from  this  boundary  should  be  construed  as  originating  from  the 
region  immediately  downstream  from  this  boundary.  The  top  and  downstream 
I  boundaries  are  simple  transmissive  boundaries  designed  to  make  a  zero-gradient 

match  to  local  flow  conditions  for  pressure  and  normal  velocity.  This  type  of 
boundary  Invariably  generates  some  sort  of  artificial  wave  reflection  upon  the 
arrival  of  waves  from  the  interior  of  the  computational  grid.  It  has  no 
resetting  or  self-corrertlng  capability,  so  there  is  no  reason  to  expect  that 
it  will  return  to  some  "correct"  value  after  being  perturbed,  even  if  the 
i  computation  is  carried  out  for  a  long  time  relative  to  the  hydrodynamic  events 

being  simulated. 

B.  Transmissive  Boundary  Check. 

The  computations  labelled  SS/TBl  through  SS/TB5  were  used  to  quantify  the 
effects  of  moving  only  the  top  transmissive  boundary  in  the  2-D  Cartesian 
grid.  ("TB"  refers  to  Top  Boundary  and  "n"  to  the  computation  number.)  The 
location  of  the  upstream  constant-inflow  boundary  was  kept  constant  at  24.37  m 
from  the  front  face  of  the  target.  The  location  of  the  downstream  simple 
transmissive  boundary  was  kept  constant  at  22.98  m  from  the  back  face  of  the 
target.  SS/TBl  was  run  for  610  ms  simulated  real  time,  with  the  top 
transmissive  boundary  of  the  grid  located  12.51  m  from  the  top  flat  section  of 
the  target.  SS/TB2  was  run  for  300.2  ms  (top  boundary  at  27.03  m) ,  SS/TB3  for 


300.0  ms  (top  boundary  at  41.55  m) ,  SS/TB4  for  610.1  ms  (top  boundary  at  56.08 
m) ,  and  SS/TB5  was  run  for  610.2  ms  (top  boundary  at  85.12  m).  SS/EB5  was  run 
for  610.1  ms,  with  a  top  boundary  location  equal  to  that  of  SS/TB5,  but  with 
the  upstream  inflow  boundary  located  135.30  ra  from  the  front  face  and  the 
downstream  boundary  located  128.63  m  from  the  back  face.  These  data  and  other 
descriptive  data  are  summarized  in  Table  1. 


The  computational  flow  field  cells  which  comprise  target  SS  and  its 
immediate  vicinity  range  from  26.40  to  22.47  cm  in  the  direction  of  shock 
travel  and  are  15.24  cm  in  height.  Overall,  the  target  was  resolved  by  12 
cells  in  depth  and  12  cells  in  height.  Each  succeeding  computation  in  the  TBn 
series  described  above  contained  the  previous  computational  flow  field  as  a 
subset.  SS/TBl  had  a  total  of  86  flow  field  cells  in  the  direction  of  shock 
travel  and  48  cells  in  height,  SS/TB2  contained  86  by  68,  SS/TB3  86  by  88, 
SS/TB4  86  by  108,  SS/TB5  86  by  148,  and  SS/EB5  246  by  148.  SS/EB5  contained 
SS/TB5  as  a  subset. 

One  of  the  first  considerations  in  the  design  of  a  grid  is  to  determine 
the  probable  amount  of  disturbance-free  time  in  the  region  of  interest  that 
can  be  simulated.  After  the  planar  incident  shock  strikes  the  front  face  of 
the  target,  a  strong,  planar,  reflected  shock  Is  sent  upstream  at  a  velocity, 
V^,  relative  to  the  oncoming  flow.  The  gas  immediately  behind  the  reflected 

shock  just  after  its  creation  is  at  rest.  An  expansion  or  weak  compression 
wave  in  this  region  travels  at  the  local  speed  of  sound,  c^.  If  the  extreme 

case  is  considered  where  there  is  no  attenuation  of  this  reflected  shock, 
hereinafter  called  Upstream  Assumption  1  (UAl) ,  then  the  travel  time,  t^^j, 

required  to  go  from  the  front  face  to  the  upstream  boundary  is, 

tufl=^- 

''r  ''s 

Here,  d^  is  the  distance  from  the  front  face  to  the  upstream  boundary  and  v^ 

Is  the  particle  velocity  of  the  gas  that  has  been  processed  by  the  incident 
shock.  The  travel  time,  required  to  go  downstream  from  the  upstream 

boundary  to  the  front  face  is 

d  f 

tdfl  =  ~  •  (2) 

Assumption  UAl  predicts  an  upstream  travel  time  of  73  ms  and  a  downstream 
travel  time  of  49  ms  for  a  total  disturbance-free  time  of  122  ms  for  the  TBn 
series.  The  comparable  times  for  SS/EB5  are  405  ms  upstream  and  272  ms 
downstream  for  a  total  travel  time  of  677  ms. 

The  more  reasonable  assumption  is  that  the  reflected  shock  will  decay  in 
strength  as  it  expands  upstream  from  the  target.  Upstream  Assumption  2  (UA2) 
adopts  the  extreme  position  that  the  effect  of  the  reflected  shock  on  wave 
propagation  between  the  upstream  boundary  and  the  front  face  is  negligible. 
Therefore,  the  travel  time,  (f*"ont  face  to  upstream  boundary)  is 
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Table  1.  Hydrocode  Computations 
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where  is  the  local  speed  of  sound  behind  the  incident  shock.  The  travel 
time,  (upstream  boundary  to  front  face)  is 

_  df  (4) 

tdf2  -  C3+V,  • 

Assumption  UA2  predicts  an  upstream  travel  time  of  296  ms  and  a  downstream 
travel  time  of  32  ms  for  a  total  travel  time  of  328  ms  for  the  TBn  series. 

The  comparable  times  for  SS/EB5  are  1.644  s  upstream  and  0.177  s  downstream 
for  a  total  travel  time  of  1.821  s. 


Assumptions  UAl  and  UA2  represent  the  two  extremes  of  the  possible  arrival 
time  at  the  front  face  of  a  disturbance  travelling  this  out-and-return  path. 

It  seems  logical  that  the  actual  disturbance-free  time  must  lie  between  these 
two  extremes.  The  two  assumptions  may  be  combined  into  a  more  general 
assumption,  UA3,  that  the  reflected  shock  makes  some  fractional  modification 
to  the  average  wave  propagation  velocities  in  the  region  upstream  from  the 
front  face.  It  is  proposed  that  the  travel  time,  (front  face  to 

upstream  boundary)  is 

_ £f _ _ 

tuf3  "  xVr+(l-x)cs-Vs  ’  (5) 

and  the  travel  time,  (upstream  boundary  to  front  face)  is 

_ ^f 

tdf3  "  xcr+(l-x)(cs+vs) 

The  factor,  x,  in  Equations  (5)  and  (6)  represents  the  fractional  contribution 
of  the  reflected  shock  to  the  upstream  wave  velocities  and  is  limited  to  a 
range  of  values  from  zero  to  one.  The  limiting  case  where  x  =  1  is  consistent 
with  UAl  and  x  =  0  is  consistent  with  UA2.  The  value  of  x  will  depend  upon 
the  geometry  of  both  the  target  and  the  computational  grid. 


The  reflected  shock  in  this  2-D  problem  very  quickly  Cakes  on  the  shape  of 
a  right-circular  cylinder.  The  surface  of  a  cylinder  varies  linearly  with 
radius,  and  the  travel  from  the  target  face  to  the  upstream  boundary 
represents  one  half  of  the  total  distance  travelled.  The  initial  "radius,"  r, 
of  the  reflected  shock  is  taken  as  the  radius  of  a  semi-circle  having  the  same 
area  as  the  cross-sectional  area  of  the  target.  The  final  radius  of  the 
reflected  shock  is  taken  as  (d^  +  r).  The  Initial  radius  divided  by  the  final 

radius  is  taken  as  the  contribution  of  the  reflected  shock,  x,  where 


X 


r 

df+r 


(7) 


Equation  (5)  yields  an  upstream  travel  time  of  254  ms  and  Equation  (6)  yields 
a  downstream  travel  time  of  34  ms  for  a  total  predicted  time  of  278  ms  for 
disturbances  to  travel  from  the  front  face  to  the  upstream  boundary  and 
return.  This  applies  to  all  computations  in  the  TBn  series.  The  comparable 
times  for  SS/EB5  are  1.581  s  upstream  and  178  ms  downstream  for  a  total  travel 
time  of  1.759  s. 
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The  disturbance-free  time  fur  the  back  face  (relative  to  artificial  wave 
arrivals  from  the  downstream  boundary)  may  be  computed  more  easily  than  that 
for  the  front  face.  The  travel  time,  t.,  ,  of  the  incident  shock  from  the 

uD 

plane  of  the  front  face  (t  -  0  there)  to  the  back  boundary  is 


db+^^tgt 

tdb  -  ‘ 

where  d,  is  the  distance  from  the  back  face  to  the  downstream  boundary,  V 
D  S 

the  incident  shock  velocity,  and  is  the  target  depth.  The  travel  time 

t  ,  ,  (back  boundary  to  back  face)  is 

UD 


(8) 

is 


<^b 

•"ub  =  cg-vs  ‘ 


(9) 


Equation  (8)  yields  a  downstream  travel  time  of  43  ms  and  Equation  (9)  yields 
an  upstream  travel  time  of  279  ms  for  a  total  disturbance-free  time  of  322  ms 
for  the  back  face  relative  to  the  downstream  boundary.  These  times  apply  to 
all  computations  in  the  TBn  series.  The  comparable  times  for  SS/EB5  are  217 
ms  downstream  and  1.563  s  upstream  for  a  total  disturbance-free  time  of  1.780 
s. 


An  analysis  similar  to  that  for  the  front  face  can  be  made  for  predicting 
the  arrival  time  at  the  top  face  of  the  target  of  a  disturbance  from  the  top 
boundary.  A  lower  limit  can  be  established  for  the  travel  time  of  dis¬ 
turbances  between  the  top  face  of  the  target  (taken  for  simplicity  to  be  the 
center  of  the  top  flat  section)  and  the  top  boundary.  Top  Assumption  1  (TAl) 
is  that  the  time,  t  j,  (top  face  to  top  boundary  and  return)  is 


where  d^  is  the  distance  from  the  top  face  to  the  top  boundary.  Equation  (10) 

predicts  a  lower-limit  disturbance-free  time  of  59  ms  for  SS/TBl ,  128  ms  for 
SS/TB2,  196  ms  for  SS/TB3,  265  ms  for  SS/TB4,  and  402  ms  for  SS/TB5  and 
SS/EB5. 

Making  the  assumption  that  there  is  no  contribution  from  the  reflected 

shock,  waves  moving  between  the  top  boundary  and  top  face  will  have  a  velocity 

component  equal  to  v  imposed  on  them,  thereby  Increasing  their  required 

s 

travel  time  and  distance.  Top  Assumption  2  (TA2)  is  that  the  disturbances 
likely  to  arrive  first  at  the  top  face  (when  considered  alone)  are  those  which 
have  originated  at  the  top  face  and  have  travelled  upstream  at  an  angle 
(relative  to  the  normal  vector  of  the  top  flat  section)  whose  tangent  may  be 
computed  by  ^ 

The  time,  required  for  a  disturbance  to  travel  to  this  limiting  upstream 

point  on  the  top  boundary  is 


dt 


2 


(df )tan9 


( 1  7) 


(13) 


The  time , 


dt2’ 


required  for  a  disturbance  to  return  from  this  point 


dt 

tdt2  =  7" 
^8 


dttan6 


is 


Equations  (11),  (12),  and  (13)  yield  combined  travel  times  of  195  ms  for 
SS/TBl,  421  ms  for  SS/TB2,  647  ms  for  SS/TB3,  873  ms  for  SS/TB4,  and  1.326  s 
for  SS/TB5  and  SS/EB5. 


A  more  probable  assumption  would  be  similar  to  that  used  for  the  upstream 
timing  analysis.  Top  Assumption  3  (TA3)  is  that  there  is  a  fractional 
contribution,  y,  from  the  reflected  shock,  where 


y  = 


df+r 


(14) 


Equations  (11),  (12),  and  (13)  are  modified  as  follows,  with  all  terms 
retaining  their  previous  definitions  and  y  limited  to  a  range  of  values  from 
zero  to  one. 


(tan9)3  = 


,(l-y) 


Cs(l-y)+yCr 


(15) 


dt 

Cg(l-y)+ycr 


+ 


d^tanB 


(16) 


tdt3 


_ it_ _ 

Cs(l-y)+ycr 


+ 


dj.tan0 

^Cs+Vs)(l-y)+yCi. 


(17) 


A  value  of  y  =  0  is  consistent  with  TA2,  and  a  value  of  y  =  I  is  the 
degenerate  case  where  a  shock  strikes  the  closed  end  of  a  constant-area  tube. 
For  SS/TBl,  y  =  0.124  and  the  angle  from  Equation  (15)  is  34.63  degrees 
measured  counterclockwise  from  the  vertical.  Equations  (15),  (16),  and  (17) 
yield  combined  travel  times  of  134  ms  for  SS/TBl,  341  ms  for  SS/TB2,  558  ms 
for  SS/TB3,  779  ms  for  SS/TB4,  and  1.226  s  for  SS/TB5  and  SS/EB5.  The  arrival 
times  predicted  by  the  various  methods  of  the  artificial  waves  from  the 
transmissive  boundaries  are  shown  in  Table  2. 

Figure  2  shows  the  average  overpressure  on  the  front  face  of  target  SS  for 
SS/TBl  through  SS/TB5.  Also  shown  in  Figure  2  is  a  curve  for  SS/EB5.  SS/EB5 
will  be  used  as  a  reference  computation  to  assist  in  identifying  the  arrivals 
of  waves  from  the  top  boundary  in  the  TBn  series  of  computations,  and  will 
Itself  later  be  used  to  illustrate  a  significant  defect  in  HULL'S  ability,  at 
least  for  a  non-decaying  shock  and  with  this  flow-field  resolution,  to 
properly  model  the  strong  vortex  to  be  expected  behind  a  blunt  2-D  target  such 
as  this  one.  (The  late  time  loading  in  SS/EB5  is  too  high  for  all  faces 
except  the  front  face  and  possibly  the  front  slope  because  of  the  improper 
modelling  of  this  vortex,  which  causes  an  excessively  strong  flow  to  impinge 
on  the  back  face  of  SS.) 


Computation  SS/TBl  shows  a  divergence  from  both  SS/TBS  and  SS/EB5  caused 
by  a  compression  wave  arrival  from  the  top  boundary  at  110  ms.  (A  wave 
arrival  Is  considered  as  confirmed  once  Che  average  overpressure  values  from  a 
computation  have  deviated  by  0.5  percent  or  more  from  those  for  a  reference 
computation.)  Until  that  time,  all  computations  are  In  agreement.  Compu¬ 
tations  SS/TB2,  SS/TB3,  and  SS/TBA  all  show  agreement  with  SS/TB5  throughout 
their  respective  run  times.  The  next  deviations  from  computation  SS/EB5  are 
by  SS/TB2  at  270  ms  and  by  SS/TB3,  SS/TB4,  and  SS/TB5  at  290  ms,  due  primarily 
to  compression  wave  arrivals  from  the  upstream  boundary.  Except  for  an 
obvious  Interference  from  Che  top  boundary  for  SS/TBl,  Che  artificial  wave 
arrivals  correlate  well  with  Che  predicted  time  of  278  ms  from  UA3  for  the  TBn 
series.  It  appears  that  SS/EBS  provides  a  reliable  estimate  for  both  the 
diffraction  phase  and  drag  phase  average  overpressure  on  the  front  face. 

Table  3  provides  a  summary  of  Che  first  deviations  of  SS/TBl  through  SS/TBA 
from  SS/TB5  for  this  face  and  all  other  faces.  These  deviations  are  direct 
Indicators  of  the  effects  of  the  locations  of  the  Cop  boundary.  Table  A 
provides  a  similar  summary  of  the  first  deviations  of  Che  TBn  series  from 
SS/EB5  for  all  faces,  thus  showing  Che  superimposed  effects  of  the  locations 
of  the  upstream  and  downstream  boundaries. 

Table  3.  First  Deviation  from  SS/TB5 


Front 

Front 

Top  Flat 

Back 

Back 

Computation 

Face 

Slope 

Section 

Slope 

Face 

(ms) 

(ms) 

(ms) 

(ms) 

(ms) 

SS/TBl 

110 

100 

80 

75 

75 

SS/TB2 

* 

270 

180 

180 

170 

SS/TB3 

* 

* 

270 

290 

290 

SS/TBA 

* 

A90 

390 

390 

390 

*Agreement 

through  Its  end  of 

run  time. 

Table  A.  First 

Deviation 

from  SS/EB5 

1  Reference) 

Front 

Front 

Top  Flat 

Back 

Back 

Computation 

Face 

Slope 

Section 

Slope 

Face 

(ms) 

(ms) 

(ms) 

(ms) 

(ms) 

SS/TBl 

110 

100 

80 

75 

75 

SS/TB2 

270 

290 

130 

130 

120 

SS/TB3 

290 

250 

130 

130 

120 

SS/TBA 

290 

250 

130 

130 

120 

SS/TB5 

290 

250 

130 

130 

120 

Figure  3  shows  the  average  overpressure  on  the  front  slope.  There  Is  more 
variation  evident  In  these  curves  because  of  the  orientation  of  the  front 
slope  relative  to  the  top  boundary.  SS/TBl  shows  an  initial  slight  deviation 
from  SS/TB5  at  100  ms,  a  return  Co  relative  agreement,  then  a  permanent 
deviation  starting  at  200  ms.  The  deviations  at  100  and  200  ms  appear  to  be 
due  to  the  arrival  of  compression  waves  from  the  top  boundary.  There  Is  also 
an  arrival  of  an  artificial  expansion  wave  at  AAO  ms.  SS/TB2  deviates  from 
SS/TB5  at  270  ms,  not  long  before  its  ending  of  300.2  ms,  and  SS/TB3  matches 
SS/TB5  throughout  Its  300  ms  run  time.  SS/TBA  does  not  diverge  from  SS/TB5 
until  A90  ms.  This  can  only  be  due  to  a  difference  in  top  boundary  effects  as 
far  as  these  two  computations  are  concerned.  However,  there  Is  more  here  chat 
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Figure  3*  Boundary  Check  Computations  for  Average 
Overpressure  on  the  Front  Slope. 


la  of  Interest.  The  only  differences  between  SS/TB5  and  SS/EB5  are  the 
locations  of  the  upstream  and  downstream  boundaries.  The  Joint  divergence  of 
SS/TB4  and  SS/TBS  from  SS/EB5  for  the  front  slope  occurs  at  250  ms,  so  this 
divergence  must  therefore  be  due  to  an  artificial  compression  wave  from  the 
upstream  boundary.  SS/EB5  also  shows  a  distinct  wave  arrival  at  290  ms,  with 
a  small  cusp-like  feature  at  310  ms.  From  the  earlier  discussion  on  the 
timing  of  disturbance  arrivals,  it  Is  very  unlikely  to  have  originated  from 
any  of  the  transmissive  boundaries.  It  is  possible  that  It  Is  from  a 
recompression  following  an  over-relief  of  the  initial  shock  reflection,  but 
that  is  speculative,  and  the  timing  seems  late  for  that.  At  present,  the  best 
that  can  be  said  is  that  It  probably  originates  from  an  unknown  upstream 
source. 

Figure  A  shows  the  set  of  curves  for  average  overpressure  on  the  top  flat 
section.  SS/TBl  shows  the  arrival  of  an  expansion  wave  from  the  top  boundary 
at  80  ms,  followed  by  some  combination  of  compression  waves  arriving  at  about 
160  ms  from  the  top  boundary  upstream  from  the  target.  These  compression 
waves  plus  others  from  the  upstream  and  downstream  boundaries  form  an 
artificial  peak  of  159  kPa  In  the  overpressure  curve  by  410  ms,  27  percent 
above  the  125  kPa  overpressure  at  the  same  time  for  SS/EB3.  SS/TB2  shows  the 
arrival  of  a  weak  expansion  wave  from  the  top  boundary  at  180  ms,  fortuitously 
cancelled  by  an  artificial  compression  wave  so  that  it  remained  close  to  the 
curve  for  SS/EB5.  SS/TB3  shows  only  the  beginning  of  an  artificial  expansion 
'.ave  at  270  ms,  and  SS/TB4  shows  a  weak  expansion  wave  arrival  at  390  ms. 

Thus,  there  Is  a  definite  correlation  between  the  expansion  wave  arrival  in 
the  TBn  series  and  the  top  transmissive  boundary  location.  A  comparison 
between  SS/TB5  and  SS/EB5  shows  that  a  weak  compression  wave  influences  the 
results  for  SS/TB5  between  130  ms  and  320  ms.  This  is  also  true  for  both 
SS/TB3  and  SS/TB4.  The  strong  compression  wave  at  320  ms  is  clearly  due  to  a 
reflection  from  the  downstream  boundary.  However,  the  compression  wave 
arrival  at  130  ms  was  initially  quite  puzzling  because  the  preceding  analysis 
showed  no  realistic  probability  of  wave  arrivals  from  the  transmissive 
boundaries  at  such  an  early  time.  The  key  to  understanding  this  compression 
wave  arrival  is  in  the  artificial,  unbounded  downstream  growth  of  the  vortex 
at  the  back  face.  (This  will  be  discussed  in  detail  later  In  the  section  on 
SS/EB5.  Only  the  specifics  necessary  to  understand  this  wave  arrival  will  be 
discussed  here.)  The  artificially  large  back  face  vortex  lengthens  In  the 
downstream  direction  at  120  m/s  (the  rate  of  progression  of  Its  downstream  end 
with  time),  while  Its  upstream  end  remains  at  the  back  face.  It  generates  an 
excessively  strong  reverse  flow  along  the  ground  plane,  having  a  velocity  of 
200  m/s  directed  at  the  back  face.  A  relatively  straightforward  timing  check 
indicates  that  the  incident  shock  arrives  at  *:he  downstream  boundary  at  43  ms, 
sends  an  artificial  compression  wave  upstream  which  meets  the  end  of  the 
vortex  at  109  ms  at  a  point  12  m  downstream  from  the  back  face,  and  rides 
upstream  on  the  reverse  flow  of  the  vortex,  arriving  at  the  back  face  at  128 
ms.  This  Is  essentially  an  exact  match  in  time  with  the  observed  disturbance 
at  130  ms  on  the  top  flat  section.  This  disturbance  was  not  seen  In  SS/TBl 
because  waves  from  the  top  boundary  washed  out  that  effect;  It  is  missing  from 
SS/EB5  because  of  the  remote  downstream  boundary.  There  is  also  a  slight 
cusp-like  feature  in  SS/EB5  at  310  ms,  essentially  at  the  same  time  as  the 
cusp  for  the  front  slope  average  overpressure  for  SS/EB5  shown  in  Figure  3, 
but  noticeably  weaker.  This  lends  support  to  the  conclusion  that  the  wave 
which  produced  it  originates  from  an  upstream  source. 
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Figure  4.  Boundary  Check  Computations  for  Average 
Overpressure  on  the  Top  Flat  Section. 


Figure  5  shows  the  average  overpressure  versus  time  on  the  back  slope  of 
target  SS«  All  computations  show  an  Initial  pressure  peak  of  112.9  kPa  at  6 
ms,  45.4  percent  below  the  Incident  shock  overpressure  of  206.8  kPa  because  of 
the  30  degree  negative  angle  of  the  back  slope  relative  to  the  undisturbed 
flow.  This  peak  Is  followed  by  a  sharp  drop  to  a  minimum  value  of  27.2  kPa  at 
15  ms  because  of  the  vortex  that  forms  off  the  trailing  edge  of  the  flat 
section  on  the  top  face.  All  computations  are  In  agreement  through  80  ms, 
when  SS/TBl  shows  the  arrival  of  the  artificial  expansion  wave  from  the  top 
boundary.  The  subsequent  performance  of  SS/TBl  is  similar  to  that  on  the  top 
flat  section,  as  are  the  rest  of  the  curves  in  the  TBn  series.  The  other 
Initial  artificial  expansion  wave  arrivals  are  at  180  ms  for  SS/TB2,  290  ms 
for  SS/TB3,  and  390  ms  for  SS/TB4.  The  comments  concerning  the  compression 
wave  arrivals  at  130  ms  and  320  ms  that  were  made  for  the  top  flat  section 
also  apply  here.  It  Is  also  interesting  that  SS/EB3  for  the  back  slope  does 
not  show  the  cusp-llke  feature  at  310  ms,  shown  clearly  for  the  front  slope 
and  weakly  for  the  top  flat  section.  This  confirms  its  upstream  origin.  The 
artificial  wave  arrivals  relative  to  SS/TB5  shown  for  the  front  slope,  top 
flat  section,  and  back  slope  indicate  that  TAl  gives  too  low  an  estimate  of 
wave  arrival  time  and  TA3  too  high  an  estimate,  with  TAl  somewhat  closer  to 
the  actual  values  on  average. 

Figure  6  shows  the  average  overpressure  on  the  back  face  of  the  target. 

All  computations  agree  through  the  first  75  ms,  showing  a  relatively  sharp 
rise  to  a  local  maximum  of  101.1  kPa  at  20  ms,  and  then  a  brief  drop  to  a 
local  minimum  ot  98.8  kPa  at  30  ms.  SS/TBl  shows  the  arrival  of  an  expansion 
wave  at  75  ms,  completing  an  interesting  progression  of  initial  wave  arrivals 
on  the  various  faces  of  target  SS:  expansion  waves  at  75  ms  for  the  back  face, 
75  ms  for  the  back  slope,  and  80  ms  for  the  top  flat  section;  and  compression 
waves  at  100  ms  for  the  front  slope  and  110  ms  for  the  front  face.  Thus,  the 
initial  disturbances  due  to  compression  waves  on  the  front  face  and  front 
slope  originate  from  a  region  of  the  top  boundary  upstream  from  the  target. 

The  Initial  disturbances  due  to  expansion  waves  on  the  top  flat  section,  the 
back  slope,  and  the  back  face  originate  from  a  region  of  the  top  boundary 
closer  to  the  target. 

The  initial  expansion  wave  which  arrives  at  the  back  face  of  target  SS  in 
SS/TBl  at  75  ms  is  followed  by  a  recompresslon  starting  at  150  ms.  This 
recorapresslon  peaks  at  201  kPa  overpressure  by  410  ms  before  it  drops  due  to 
another  expansion  wave,  most  likely  originating  from  the  top  boundary.  (This 
strong  expansion  wave  is  absent  in  SS/TB4  and  SS/TB5.)  The  recompresslon 
before  350  ms  is  due  to  compression  waves  from  the  top  boundary  and  the 
downstream  vortex  action  discussed  earlier,  and  after  350  ms  is  due  to  a 
compression  wave  from  the  Incident  shock  interaction  with  the  downstream 
boundary.  SS/TB2  for  the  back  face  of  target  SS  in  Figure  6  shows  a  weak 
expansion  wave  arrival  at  170  ms,  10  ms  ahead  of  the  arrival  of  the  Initial 
expansion  wave  for  the  back  slope,  which  is  consistent  with  the  pattern  noted 
for  SS/TBl.  SS/TB3  shows  a  weak  initial  expansion  wave  at  290  ms,  and  SS/TB4 
shows  an  expansion  wave  arrival  relative  to  SS/TB5  at  390  ms.  Both  compu¬ 
tations  show  the  arrival  of  a  very  weak  expansion  wave  at  120  ms  relative  to 
SS/EB5.  The  flow  on  the  back  face  has  become  exaggerated  by  this  time,  as 
will  be  shown  later,  so  this  may  not  be  outside  of  the  realm  of  possibility. 
The  flow  is  recompressed  by  the  wave  moving  up  the  vortex  starting  at  160  ms 
and  by  a  strong  compression  wave  from  the  downstream  boundary  at  350  ms.  In 
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the  absence  of  top  boundary  effects,  the  predicted  322  ms  arrival  time  of 
artificial  waves  from  the  downstream  boundary  for  the  TBn  series  appears  to  be 
of  acceptable  accuracy. 

C.  Extended  Boundary  Reference  Computation  SS/EB5. 

As  indicated  earlier,  although  SS/EB5  is  useful  as  a  reference  computation 
for  the  TBn  series  in  the  previous  subsection,  its  modelling  of  the  vortex 
behind  target  SS  at  late  time  is  not  correct.  Figure  6,  for  example,  shows 
that  the  average  overpressure  on  the  back  face  from  SS/EB5  is  steadily 
increasing  with  time.  The  ending  time  of  610.1  ms  for  SS/EB5  represents  121.5 
shock  crossing  times.  (A  shock  crossing  time  is  defined  as  the  time  required 
for  the  incident  shock  to  travel  from  the  plane  of  the  front  face  to  the  plane 
of  the  back  face.)  This  number  of  shock  crossing  times  should  be  more  than 
sufficient  for  a  steady  state  loading  to  be  reached,  yet  SS/EB5  still  shows  a 
significant  positive  slope  in  the  average  overpressure  versus  time  curve  at 
its  ending  time. 

Figure  7a  shows  a  close-up  view  of  the  velocity  vectors  at  10  ms  in  the 
flow  field  in  the  vicinity  of  target  SS.  This  is  relatively  soon  after  the 
shock  has  passed  over  the  target,  and  the  downstream  recirculation  region  has 
not  yet  developed.  Figure  7b  shows  the  same  region  at  20  ms.  The  recir¬ 
culation  region  has  just  begun  to  form,  with  its  center  of  rotation  located 
0.7  m  downstream  from  the  back  face.  Its  scale  and  location  appear  to  be 
quite  reasonable.  Figure  7c  shows  the  downstream  recirculation  region  fully 
formed  by  30  ms,  with  its  center  1.7  m,  or  0.56  target  depths,  downstream  from 
the  front  face.  The  reverse  flow  near  the  ground  plane  is  beginning  to 
impinge  strongly  on  the  back  face,  raising  the  pressure  there  and  also  on  the 
back  slope.  The  pressure  rise  on  the  back  slope  is  high  enough  that  the  flow 
is  turned  away  from  it  and  is  nearly  parallel  to  the  ground  plane.  By  40  ms. 
Figure  7d,  the  center  of  the  recirculation  region  has  moved  to  2.6  m  down¬ 
stream  from  the  back  face.  The  recirculation  region  has  an  overall  length  of 
4.5  m.  The  upstream  counterflow  is  now  impinging  on  the  back  face  with 
sufficient  strength  that  the  pressure  on  the  back  face  and  back  slope  is  high 
enough  that  the  flow  above  the  back  slope  has  a  significant  positive  vertical 
component.  This  seems  to  be  the  beginning  of  the  artificially  large  down¬ 
stream  vortex.  The  vortex  is  centered  at  2.8  ra  with  an  overall  length  of  5.0 
m  by  43  ms  (Figure  7e),  centered  at  2.9  ra  with  a  length  of  5.5  m  by  46  ms 
(Figure  7f),  and  centered  at  3.4  m  with  a  length  of  about  6.0  m  by  50  ms 
(Figure  7g).  It  is  during  this  time  from  40  to  50  ms  that  the  computation  in 
this  region  of  the  flow  field  loses  credibility.  The  downstream  vortex  seems 
to  grow  without  bound  in  the  downstream  direction,  with  a  constant  height  of 
approximately  2.5  ra.  Figure  7h  shows  the  local  region  of  the  flow  field  at 
the  ending  time  of  610.1  ms.  The  "vortex"  has  evolved  into  a  high-speed 
upstream  flow  along  the  ground  plane,  a  low-speed  layer  above  it,  and  then  a 
top  layer  blending  into  the  downstream  flow  above  it,  all  extending  tens  of 
metres  downstream  from  the  target. 

It  seems  unlikely  that  the  center  of  the  downstream  recirculation  region 
at  the  back  face  would  extend  more  than  the  maximum  dimension  of  the  target, 
3.05  m  in  depth  for  target  SS.  If  this  is  accepted  as  a  plausible  criterion, 
then  the  latest  time  for  "acceptable"  results  for  the  back  face  and  back  slope 
is  46  ms  (see  Figure  7f).  Figure  8  shows  the  drag  coefficient  computed  for 
target  SS  using  the  uncorrected  loading  for  the  front  and  back  faces  and 
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slopes  from  SS/EB5,  and  the  "corrected"  drag  coefficient.  The  corrected  drag 
coefficient  was  computed  by  assuming  that  the  front  face  and  front  slope 
loadings  are  acceptable  up  to  290  ms.  and  then  constant  at  the  290  ms  value 
for  the  remainder  of  the  computation  time  to  eliminate  the  effect  of  the 
artificial  wave  arrival  discussed  earlier.  The  loadings  on  the  back  face  and 
back  slope  were  assumed  to  be  acceptable  through  46  ms,  then  constant 
thereafter. 

The  drag  coefficient  as  It  is  used  here  is  actually  an  area-weighted 
pressure  coefficient,  taking  the  loading  due  to  overpressure  on  the  front  as 
positive  and  that  on  the  back  as  negative.  This  Is  also  the  sense  in  which 
the  presentation  of  the  drag  coefficient  for  the  diffraction  phase  as  shown  in 
Figure  8  should  be  taken.  In  this  way,  a  curve  that  is  complete  In  time  is 

presented.  This  was  also  done  a  number  of  years  ago  by  Tucker^  in  an  eval- 
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nation  of  an  earlier  study  describing  experimental  drag  coefficients  for 
bluff  bodies  struck  by  non-decaying  shocks. 

The  curve  for  drag  coefficient  in  Figure  8  using  the  uncorrected 
overpressure  falls  steadily  from  its  peak  value  of  4.70  during  the  diffraction 
phase,  reaching  the  unacceptably  low  value  of  0.96  at  the  610.1  ms  ending 
time.  The  curve  using  the  corrected  values  matches  the  uncorrected  curve 
through  46  ms;  a  steady,  corrected  drag  coefficient  of  1.42  is  reached  at  290 
ms.  This  value  seems  to  be  plausible  for  this  target,  especially  considering 
that  HULL  does  not  model  viscosity  at  all.  Some  drag  coefficients  for  simple 
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shapes  published  by  Hoerner  may  be  used  as  comparative  guidelines.  A  plate 
mounted  between  parallel  walls  and  oriented  normal  to  a  Mach  0.8  flow  (the 
local  Mach  number  behind  the  shock  in  the  present  study)  has  a  drag  co¬ 
efficient  of  2.2  (Reference  9,  p. 16-16),  a  2-D  square  cross-section  cylinder 
on  a  ground  plane  has  a  drag  coefficient  of  2.05  (ibid.,  p.3-l7),  and  a 
right-circular  cylinder  In  a  Mach  0.8  cross-flow  has  a  drag  coefficient  of  1.5 
(ibid.,  p. 16-16).  The  plate  and  square  cylinder  represent  upper  limits  for 
the  drag  coefficient  for  target  SS,  and  the  circular  cylinder  represents  a 
lower  limit.  If,  for  example,  the  drag  coefficient  computed  from  the 
corrected  values  for  SS/EB5  is  assumed  to  be  too  low  by  10  to  20  percent 
(because  of  either  undetected  or  uncorrected  errors  or  simply  the  lack  of 
viscous  effects),  then  a  quite  consistent  drag  coefficient  of  from  1.56  to 
1.70  is  implied.  Thus,  it  seems  possible  to  estimate  drag  loading  on  targets 
with  HULL,  but  only  if  great  care  Is  taken  in  Interpreting  the  results.  The 
far  more  preferable  choice  would  be  to  use  a  hydrocode  that  is  more  complete 
in  its  approximation  of  the  Navier-Stokes  equations,  particularly  in  regard  to 
the  Inclusion  of  cross-derivative  terms  for  better  handling  of  rotational 
flows,  and  the  inclusion  of  viscous  effects. 

D.  Cross-Check  of  Computations. 

The  next  set  of  figures  shows  comparisons  between  SS/TBl  and  three  similar 
computations:  SS/3-D  (Stair-Stepped  target,  3-D  coding),  RT/L-R  (Ramped 
Target,  Low  Resolution),  and  RT/H-R  (Ramped  Target,  High  Resolution).  SS/3-D 
uses  the  same  cross-section  of  target  SS  as  used  for  the  TBn  series  and  is 
made  hydrodynamically  2-D  by  using  the  3-D  coding  with  a  five  cell  wide  grid 
having  perfectly-reflecting  end  planes;  it  has  the  same  resolution  as  SS/TBl. 
Thus,  a  comparison  between  SS/TBl  and  SS/3-D  is  a  direct  Indicator  of  any 
systematic  differences  in  results  obtained  by  using  the  2-D  coding  as  opposed 
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CO  using  the  3-D  coding.  RT/L-R  uses  target  RT,  which  Is  Identical  to  target 
SS  except  that  Che  front  and  back  slopes  of  RT  are  made  smooth  through  the  use 
of  shore  cells;  It,  Coo,  has  the  same  resolution  as  SS/TBl.  A  comparison 
between  SS/TBl  and  RT/L-R  Is  a  direct  Indicator  of  the  effect  of  using  stepped 
versus  smooth  slopes.  The  grid  for  RT/H-R  was  made  by  halving  each  flow-field 
cell  In  the  grid  for  RT/L-R  In  each  direction,  a  general  doubling  of  resolu¬ 
tion  of  RT/H-R  over  RT/L-R.  Each  of  these  computations  has  problems  resulting 
from  wave  interactions  with  computational  transmissive  boundaries  which  may 
either  have  been  exacerbated,  unaffected,  or  mitigated  by  the  variations 
Introduced  In  target  generation  and  resolution;  those  aspects  of  these 
computations  will  not  be  covered  In  this  section. 

Figure  9  shows  the  average  overpressure  on  the  front  face  for  these  four 
computations.  There  Is  little  difference  between  the  computations  on  this 
flat  surface  normal  Co  Che  Incident  flow.  Figure  10  shows  Che  average  over¬ 
pressure  on  the  front  slope.  SS/TBl  and  SS/3-D  show  nearly  identical  results 
through  the  run  time  of  SS/3-D.  There  is  a  large  difference  between  RT/L-R 
and  SS/TBl,  which  Is  directly  attributable  to  modelling  the  slope  as  a  smooth 
surface  in  RT/L-R.  By  300  ms,  the  average  overpressure  on  the  front  slope  for 
SS/TBl  is  15  percent  greater  than  for  RT/L-R.  The  effects  of  shock  inter¬ 
action  with  stepped  versus  smooth  surfaces  has  also  been  documented  In  ocher 

studies.^  The  doubled  resolution  of  RT/H-R  over  RT/L-R  produces  a  lower 
average  overpressure  on  the  front  slope  for  RT/H-R.  The  average  overpressure 
on  the  top  flat  section  in  Figure  11  shows  a  pattern  similar  to  Chat  on  the 
front  slope.  RT/H-R  apparently  does  the  best  at  resolving  the  low-pressure 
vortex  that  forms  off  the  leading  edge  of  the  flat  section  at  early  time. 
Figure  12  shows  the  average  overpressure  on  the  back  slope.  Here  again  SS/TBl 
and  SS/3-D  show  close  agreement  with  each  other,  RT/L-R  and  RT/H-R  agree  very 
closely  through  the  first  20  ms  and  then  diverge,  with  RT/H-R  showing  a  higher 
value  by  130  ms.  Figure  13  shows  the  average  back  face  overpressure.  SS/TBl 
and  SS/3-D  show  very  close  agreement,  and  share  the  lowest  of  Che  three 
initial  peaks  at  101  kPa.  Interestingly,  RT/L-R  has  the  highest  initial  peak 
of  115  kPa,  and  RT/H-R  has  the  Intermediate  peak  of  110  kPa  but  has  a  value  at 
130  ms  higher  than  that  for  RT/L-R.  The  overall  effects  of  these  variations 
are  best  seen  In  Figure  14,  which  shows  the  uncorrected  drag  coefficients  for 
these  four  computations.  There  Is  essentially  no  difference  between  the  drag 
coefficients  for  SS/TBl  and  SS/3-D,  RT/L-R  tends  to  become  slightly  less  than 
SS/TBl  and  SS/3-D  at  later  time,  and  RT/H-R  is  lower  still. 

This  part  of  the  study  was  intended  to  Illustrate  the  potential  for 
differences  In  computed  loading  at  late  time  for  this  target  and  shock 
combination,  depending  on  the  method  of  target  simulation  and  flow-field 
resolution.  It  would  not  be  appropriate  to  use  these  four  computations  as  a 
direct  basis  for  applying  a  further  correction  to  the  drag  coefficient 
computed  from  SS/EB5.  The  greater  part  of  the  divergence  between  these  four 
computations  shown  In  Figure  14  occurs  after  all  four  computations  may  be 
assumed  to  have  some  combination  of  errors  due  to  either  boundary  effects  or 
Improper  modelling  of  the  downstream  vortex. 
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V.  CONCLUSION 


The  simple,  zero-gradient,  transmissive  boundaries  used  in  the  HULL 
hydrocode  cannot  be  expected  to  allow  waves  of  any  type  to  pass  through  them 
without  reflecting  an  artificial  wave  of  the  same  type  back  into  the  compu¬ 
tational  grid.  The  artificial  reflected  wave  seems  to  have  a  strength  that  is 
directly  proportional  to  the  strength  of  the  incoming  wave.  These  boundaries 
are  not  self-correcting,  because  no  solution  is  Imposed  on  them.  Thus, 
running  a  computation  for  an  extended  time  will  not  assure  the  user  that  the 
boundaries  will  return  to  some  "correct"  state.  These  artificial  disturbances 
propagate  throughout  the  computational  flow  field.  The  velocities  of  these 
disturbances  and  their  arrival  times  at  regions  of  Interest  can  be  estimated 
using  a  relatively  simple  analytical  technique  based  on  geometry,  incident 
shock  velocity,  and  the  shocked  gas  properties.  Thus  the  results  from  this 
analysis  can  be  used  to  identify  which  boundaries  or  regions  of  boundaries  are 
likely  to  be  the  source  of  particular  disturbances. 

The  analysis  of  the  computations  presented  here  was  made  more  difficult  by 
the  Improper  modelling  by  HULL  of  the  vortex  that  forms  immediately  downstream 
from  the  back  face  of  this  blunt  2-D  target  after  the  diffraction  phase.  As 
the  computation  progressed,  the  length  of  the  “vortex"  grew  seemingly  without 
bound  in  the  downstream  direction  while  maintaining  a  constant  height  roughly 
comparable  to  the  height  of  the  target.  This  generated  an  artificially  strong 
reverse  flow  along  the  ground  plane  directed  toward  the  back  face  of  the 
target.  This  caused  a  large,  additional  artificial  loading  of  the  back  of  the 
target  at  a  very  early  time,  plus  an  Increase  in  the  vertical  loading  of  the 
top  flat  section  and  back  slope.  Hence,  the  net  effect  of  the  overstated 
vortex  was  to  produce  an  artificially  increased  hold-down  force  and  an 
artificially  decreased  net  horizontal  force  as  the  drag  phase  progressed,  so 
the  resistance  of  the  target  to  translation  and  overturning  is  overstated.  It 
is  doubtful  that  the  use  of  greater  resolution  in  the  flow  field  would  do 
anything  to  make  a  fundamental  improvement  in  the  modelling  of  the  vortex.  At 
best,  the  artificial  growth  might  be  delayed,  but  probably  not  for  long.  The 
inability  to  dissipate  such  rotation  is  probably  Intrinsic  to  codes  like  HULL 
that  solve  the  inviscid  Euler  equations  for  multidimensional  flow  in  sets  of 
one-dimensional  steps  using  only  pressure  as  the  driving  potential. 

HULL  is  still  an  excellent  code  for  free  airblast  and  diffraction  phase 
loading  of  both  2-D  and  3-D  targets:  it  is  both  fast-running  and  accurate.  It 
also  does  well,  considering  that  it  is  inviscid,  for  use  in  modelling 
late-time  loading  of  simple  3-D  and  quasi-3-D  targets  because  vortex  growth  in 
the  downstream  direction  is  inhibited  by  the  transverse  pressure  gradients. 
Thus,  back  face  loading  does  not  have  the  same  opportunity  for  artificial 
growth  as  it  does  in  2-D  flow  over  a  blunt  target.  Care  must  be  taken  to  keep 
any  transmissive  boundaries  remote  from  the  region  of  interest. 
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